Mimicking general relativity in the solar system 
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In order for a modified gravity model to be a candidate for cosmological dark energy it has to 
pass stringent local gravity experiments. We find that a Brans-Dicke (BD) theory with well-defined 
second order corrections that include the Gauss-Bonnet term possess this feature. We construct the 
generic second order theory that gives, to linear order, a BD metric solution for a point-like mass 
source. We find that the Eddington parameter 7, heavily constrained by time delay experiments, 
can be arbitrarily close to the GR value of 1, with an arbitrary BD parameter uj-bu- We find the 
region where the solution is stable to small timelike perturbations. 



Brans-Dicke (BD) theory is a simple modification of 
general relativity (author?) pj] (see also (author?) ^ 
for generalisations) as it is a single massless scalar-tensor 
theory whose only parameter is the kinetic coupling term 
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Its GR limit is obtained for lJbd 00. BD gravity 
breaks the strong equivalence principle and yields at local 
scales differing Eddington parameters /3 and 7 to those 
of GR, which are strictly equal to 1. In particular the 
parameter 7, which measures how much spatial curva- 
ture is produced by unit rest mass (see (author?) (!]), 
is given by 7 = (1 -I- u;bd)/(2 -I- Wbd)- It is strongly con- 
strained by time delay experiments, such as the one con- 
ducted with the Cassini spacecraft, which recently gave 
I7 — 1| < 10~^ (author?) \^ (for a recent review and al- 
ternative methods to measure 7 see (author?) Q). This 
implies lo-b-d > 40000, therefore the scalar sector is very 
weakly coupled. 

On the other hand modification of general relativity is 
possible or even needed in order to explain effects on 
cosmological and galactic scales, or at scales just be- 
yond the solar system (such as the Pioneer anomaly (au- 
thor?) Hi). At cosmological scales, some 10^^ times big- 
ger than solar system scales, supernovae data (author?) 
[31 entertain the possibility that GR may be modified 
at large distances (author?) f^. Modification of GR is 
also being envisaged at galactic scales in order to explain 
deviations from standard Newtonian gravity in galac- 
tic rotation curves, as in MOND or Bekenstein-Sanders 
theor (author?) Scalars have been quite naturally 
introduced in order to mediate gravity modification or 
even as sources of cosmological dark energy (author?) 

These modifications are well into the classical, in- 
frared sector of gravity at very low energies, very far from 
the the Planck scale UV sector, where quantum gravity 
becomes important. Even from the point of view of UV 
modifications, string theory predicts a zoo of scalars that, 
if massless, would give BD-type phenomenology in the so- 



lar system. It is therefore quite fair to say that there is 
increasing tension between gravitational constraints im- 
posed experimentally for weak gravity in the solar system 
and laboratory tests, as well as strong gravity from bi- 
nary pulsars (see for example (author?) [lQ|), and on 
the other hand theories of modified gravity or dark en- 
ergy that aim to explain unexpected outcomes of novel 
experimental data. The aim of this letter is to show 
that well motivated second order corrections to simple 
BD theory (with no potential) can mimic GR at the so- 
lar system scale, in the sense of giving 7=1 indepen- 
dently of the parameter Wqd- This does not mean that 
the toy scalar-tensor theory in question and GR would 
not be distinguishable, quite the contrary, on cosmologi- 
cal scales their phenomenology would be totally different 
and even at the solar system level one could detect some 
effect, most probably by carrying out an experiment to 
measure (3, as we will discuss in the concluding remarks. 

Our starting point is the general (modulo field redefini- 
tions) scalar-tensor Lagrangian of second order in powers 
of the curvature tensor which has the unique property of 
giving second order field equations, 

+ 6(V(/.)2V20 + U^(t>)^ - 2V] - 16^Go£mattter • (2) 

The theory ([2]) is parametrised by the potential V and 
couplings which are all functions of the scalar field 
(j). The Gauss-Bonnet term is £gb = Rfj.iyapR'^'"^'' — 
AR^^R^^" + R^ and is a topological invariant in 4 dimen- 
sions for GR but not for scalar-tensor Theories of 
the above form ([2]) have been proposed as a solution to 
the dark energy problem (author?) [lH. For the case 
of /i = 1, it has been shown that solar system data can 
impose severe restrictions on the coupHngs (author?) 

which allows the range of possible gravity modifica- 
tion to be narrowed down. Throughout this article we 
work in units with c = 1. We will assume ^3 = ^4 = in 
the following as a compromise between introducing new 
free coupling functions in the theory and generality of 
the setup. We choose to keep the non-minimal interac- 
tion terms between the graviton and the scalar ^1 and ^2 
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rather than higher order scalar corrections. When search- 
ing for solutions we will assume also a vanishing potential 
V (as in BD) since invoking a large mass is an obvious 
but totally ad hoc way to evade solar system constraints 
and would suppress observable effects at all scales. Our 
work also contrasts with 'chameleon' models (author?) 

in which the scalar's gravity is suppressed at local 
(but not cosmological) scales by a suitably chosen poten- 
tial that yields a background dependent mass in eq. 
The higher order corrections are chosen so that the grav- 
itational propagator does not pick up extra degrees of 
freedom. Therefore the resulting field equations will be 
a-priori ghost-free around the Minkowski vacuum, second 
order in the derivatives, with well-defined Dirac distribu- 
tional terms. 

In order to derive the post-Newtonian equations, we 
now assume a static, spherically symmetric metric in 
isotropic coordinates, 



-(1 - 2U)dt^ + (1 + 2T)S^jdx'dx^ + ©(e^/^) , (3) 



where U is the Newtonian potential and T is the lead- 
ing post-Newtonian spatial contribution. They are both 
functions of the radial co-ordinate r only and are assumed 
to be of the order of the smallness dimensionless param- 
eter e = GruQ/r where uiq is the solar mass and r is a 
characteristic length scale of the problem. For the solar 
system e < 10~^ for r greater than the sun's radius. The 
post-Newtonian parameter (PPN) we will be calculating 
is Eddington's parameter defined as 7 = T/U. Matter 
energy density pm is given by the mass of the sun and is 
as usual assumed to be a distributional source at r = 0: 
Pm = TOo5(^^(a;). This is an excellent assumption given 
that the Schwarzschild radius of the sun is of the order 
of 3 km compared to scales of the order of astronomical 
units. Note however, that higher order terms in e have to 
be included in ^ in order to calculate f} (for the advance 
of Mercury's perihelion for example). For the relativistic 
experiment we will consider here, namely time delay, our 
expansion is necessary and sufficient. Let us now define 
the operators, 

i i,j 

(4) 

which will be the technical tool essential for our 
analysis. For functions with only r-dependence they 
reduce to AF = r-^dr{r^drF) and V{X,Y) = 
—2r~'^dr{rdrXdrY) and in particular. 



I?(r-",r-") = 
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and thus we can easily evaluate the relevant distribu- 
tional parts associated with V. We do not make any 



assumptions about the relative sizes of fi, ^i, V, or their 
derivatives (since we expect the higher order terms to 
play a significant role), and instead include the leading 
order in e contribution from each term in the field equa- 
tions 



fiAU 
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where we have defined ^ = d^(^2 and /2 = d^h2. We will 
also expand /i to first order in e, /i = $0 + 0{e). The 
scalar field equation on the other hand is globally of one 
order higher and gives to leading order 

90/i2A(/. + Ah2 = 2d^V -f- 2^0/1 A(2T - U) 

- Sd^^iViU, T) + ViU - T, C2) 

+ 9^C22?(C/ - T, 0) + 0{e\eV, he, e^^, • (8) 

For comparison, BD theory (V, ^ = 0) has 



/i = $ = $0 



/2 = « ^ + OicP) (9) 



Since we are expanding the equations to the lowest non- 
trivial order, we cannot estimate the second PPN param- 
eter /?, which requires higher order terms in the metric 
coefficients. This loss in generality in the metric coeffi- 
cients is compensated by the great generality of our so- 
lution. 

It is useful to define 77 = 1 — 7 whereupon the various 
constants in the model are related by Wbd = — 2 -I- l/ry 
and G = 2Go/[(2 - 7y)$o]- We see that ry = gives 
exactly GR. We are interested in gravitational theories 
emanating from ^ , which while not identical to general 
relativity, give almost identical predictions for the weak 
field of the solar system. One approach to this problem 
would be to solve the field equations of the previous sec- 
tion for a range of coupling functions a, fi, and then 
compare the resulting potentials U and T, with those of 
Einstein gravity. We will not take this approach since we 
have no interest in the solutions of the field equations, ex- 
cept for the special cases where they give (to this order in 
e) precisely the Newtonian result J7 « T « GniQ/r. This 
in particular gives us agreement with tests of Newton's 
law from planetary orbits. Therefore, instead of trying 
to find the metric ^ which solves the field equations for 
given a, fi, we will inversely start by assuming the de- 
sired Newtonian form of U and T, and then view ([G])-® 
as equations for the coupHng functions fi, ^ parametris- 
ing the theory As discussed above, we will now set 
V = 0, just as in the standard BD model and we allow 
the PPN parameter 7 = 1 — ry to take any value. Hence 
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we take 



U = 



GmQ 



T = (1 - ry) 



GmQ 



(10) 



Note also that the effective gravitational coupling G need 
not be equal to the fundamental parameter for the grav- 
itational coupling of matter Go . The general solution of 
Je])-® is then 

- 3(1 - 477 + 2Tf)S - 77(3 - 2ri)rdrS] , 
77 Ar^ /■ r'^drS 
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-dr . 



(9.0)'6 = Grr.0[2-(l-77)-^— ^9,^], (11) 
r 2 

where A is an arbitrary, dimensionful constant, obtained 
by the distributional part appearing in the equations of 
motion ([G])-® as the boundary condition at r = 0. On 
the other hand, S{r) is an arbitrary function with the reg- 
ularity condition rS' = 65* as r ^ 0, i.e. 5 = rdrS = at 
r — Q. Viewing the above expressions as the solution of 
ordinary inhomogeneous differential equations for fi, £,i, 
the function S{r) parametrises the general homogeneous 
solution of ([B])-® whereas A parametrises the particular 
solution. The integrals range from 00 to r. The gravita- 
tional coupling satisfies 



^[1-477 + 2772]. (12) 



The above equations fully specify the coupHngs needed 
to reproduce a PPN parameter 7 and an exactly New- 
tonian 1/r gravitational potential. In fact A and S now 
parametrise the theory |[2|. Setting S — X = gives us 
pure BD ^ with (p = GniQ^orj/r. Setting on top of 
that ?7 = gives GR. The key point however is that if 
we set 77 = keeping S and A non-zero we have the same 
post-Newtonian limit as standard GR, i.e. 7 = 1 in ([3]) 
without the theory actually being GR. Indeed note that 
the corresponding kinetic coupling /2 can take arbitrary 
values parametrised by A and S. Indeed when the higher 
curvature terms are included, the Newtonian potential is 
still proportional to 1/r. For non-trivial S, this is be- 
cause all the corrections to standard gravity cancel out, 
making the gravity modifications 'invisible' to this order. 
On the other hand, if A is non-zero, the corrections do not 
cancel, but instead 'mimic' Newtonian gravity. This can 
be seen from the fact that the effective gravitational cou- 
pling G receives a A dependent correction (fT2| . A similar 
effect was found for fiCcs) gravity in (author?) [l^ . 
although the resulting 7 was too large. 

This fact is made clearer when we note that the above 
solution ifTTj) does not give a specific form for (p. This is 
natural, since by a change of variables, (f) can be made 



to take any desired form. Since we wish to express the 
functions in terms of 6, let us take 



fi — ■ 

r 



(13) 



The constant (pi simply corresponds to a re-scaling of 1 
Defining = Gttiq , the expressions ([TT| then give 



/i 

6 



*o + (3'o?7 + A) 



2-77 



1 



-2[$o?? + A-2($o + A)7;2 

3(1 - 477 + 2772)5 + 77(3 - 2f])(j)d^S] . 
rlcPl 



2A,^ , , 3-277 ^ ■ 
-^(1 - 77) H :r-7^0d,S 



2(f>^ 



(14) 



= now implying 
00. S can then 
obtaining the gen- 



with the regularity conditions at r 
S and (pd^S tending to zero as (j) - 
be expanded as = J2n>i'^n'P~" , 
eral asymptotic solution to all orders in cf). Note also 
that the higher order couplings are rg dependent which 
follows from the fact that they are of dimension length 
squared. In particular this means that if we introduce the 
Gauss-Bonnet coupling constant a then it is related via 
a multiplicative number to the only length scale of the 
problem rg, namely, a = —4>ir]Xrg/8. In fact the multi- 
plicative constant is the hierarchy generated between the 
classical scale rg and y/a. 

To illustrate our result we consider the simplest case 
of S" = 0, and take I77I < 10~^ to agree with solar 
system constraints. Without loss of generality we set 
01 = A + $0??- We see that even with 77 = 

we have a BD Lagrangian ([2]) with the additional term 
£,2G^'^V ^(pV v4> which can reproduce general relativity up 
to the first post-Newtonian parameter 7. Therefore find- 
ing 7=1 does not guarantee the absence of a scalar 
interaction even in this simplest of cases since the scalar 
coupling Wbd is still freely given by <i>o/A. Note also that 
the strength of the scalar interaction Wbd is inversely pro- 
portional to the strength of the higher order corrections, 
as parametrised by A. 

We will now examine the stability of the solution lfT4|) 
with respect to time-dependent perturbations. We take 
U ^ U + SU, etc. and keep the leading order time deriva- 
tives of 6U (up to 0{SU), for linear gravity terms, and 
0{e6U) for quadratic terms). For simpHcity we will re- 
strict ourselves to the extreme case of solution (fT4|) with 
r] = and S = . The corresponding BD-like param- 
eter for the higher order theory is Wbd ~ ^o/X. The 
perturbation equations then reduce to 



3wbd<5T 



-WbdA dU ■ 



2 (pi 



(15) 
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r<pi 
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re not as robust in their GR prediction as one might 



= ASU 

01 



2rdr 



driSU-6' 



where we have assumed that the perturbations are more 
regular than the leading order solution 1/r as r 0. 
If lfT4|) is to be a viable gravity model, there needs to 
be a reasonable range of parameters for which 5U, etc. 
oscillate, rather than growing over time. Substituting 
SU{t,l^) — SU{t)expik^ and analogously for S<I),ST, 
we find 



2ujbd 



2ujBD + 3 



2cj 



BD 



9- 



18fc3r3(fcr- 
(fe^r^+4fer- 



-3) 

IF 



2w 



BD 



(18) 
the time- 



(the limit is for r — > oo). In the same limit, 
dependence of SU, ST is the same as for 6(j). 

If Wbd > 9/2 or Lo-B-D < —3/2, the perturbations will 
oscillate for large r, rather than grow exponentially, in- 
dicating that our gravitational solution mimicking GR is 
classically stable for a reasonable range of parameters, 
at least at first order. Although it might be that in- 
cluding higher order terms in the perturbation expan- 
sion the oscillations turn out to be of a growing na- 
ture, one has to be reminded that the oscillations will 
be naturally damped by the emission of gravitational 
waves and by the background cosmological expansion. 
A further constraint comes from requiring positive grav- 
itational coupling. For the above case, lfT2|) reduces to 
G = (Go/3'o)2wbd/(2u;bd - 1). Hence G > implies 
lvbd < or Wbd > 1/2, which are already covered by the 
above ranges. We expect qualitatively similar results for 
more general (fT4l) with r] or 5(0) non-zero, although a 
proof of this is beyond the scope of this paper. Let us 
remember that S{r) is anyway an arbitrary function and 
one can always choose it a posteriori in such a way to 
maintain stability. 

In this paper we exhibited a sensible second order (in 
powers of the curvature tensor not derivative) scalar- 
tensor theory which shares some characteristics of ordi- 
nary BD or GR, in particular, well-defined second order 
field equations, distributional boundary conditions and 
well defined stable vacua. We found that such a theory, 
given the right coupling functions, can mimic a GR Ed- 
dington parameter 7 exactly equal to 1 with virtually no 
constraint on the kinetic coupling ojbd (in ordinary BD 
theory actual measurements of 7 give Wbd > 40000). In 
this sense we saw that the inclusion of higher order op- 
erators in the action can mimic GR with a scalar-tensor 
theory. We do not view the solutions we have found lfT4|) 
or even the model in question |[2|, as some fundamen- 
tal scalar-tensor theory; our aim was rather to see how 
robust were the solar system predictions to higher or- 
der corrections. Our conclusion is that certain solar sys- 
tem constraints known to rule out theories such as BD 



hink. In fact similar results have been shown for cer- 
Ntain vector-tensor theories (author?) [l3| although one 
■7|xpects closer agreement with the Eddington parameters 
in the case of vectors rather than scalars. 

This does not mean that one cannot distinguish be- 
tween such higher order scalar tensor theories and GR. 
For a start the second PPN parameter (3 may not be 
unity for such theories, although if we allow for the re- 
maining higher order operators £,3,^,4 in ([2]), in principle 
we have the mathematical flexibility in the equations to 
again fix = 1 by solving for the coupling functions. The 
main difference between GR, BD and these higher order 
theories is that the coupling functions are dimensionful. 
Thus the relevant solutions such as will depend on 
the length scale of the solution, namely r^, times some di- 
mensionless number whose magnitude will determine the 
"fine-tuning" one has to impose between the length scales 
of the theory and the solar system. In other words we 
would view the experimental error bars as hierarchies be- 
tween the higher order couplings and local scales where 
the experiment is carried out. This relation may also 
be relaxed by allowing for the general second order the- 
ory at the expense of introducing further free parameters 
in the theory We further note that we have con- 
structed gravitational theories which exactly reproduce 
the Newtonian potential for the sun: U = rg/r. In fact 
it is perfectly acceptable to have U — rg/r non-zero, but 
smaller than the experimental bounds from planetary or- 
bits. The above issues as well as a calculation of /3, other 
observational signatures of such higher order theories, in 
the laboratory or in the solar system, and their cosmol- 
ogy are open interesting questions which we hope shall 
be addressed in the near future. 
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